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Abstract. Complex Hadamard matrices have received considerable attention in the past 
few years due to their appearance in quantum information theory. While a complete char- 
acterization is currently available only up to order 5 (in [5]), several new constructions of 
higher order matrices have appeared recently [H Q~21 [U [71 [TT] . In particular, the classification 
of self-adjoint complex Hadamard matrices of order 6 was completed by Beuachamp and 
Nicoara in jj], providing a previously unknown non-affine one-parameter orbit. In this paper 
we classify all dephased, symmetric complex Hadamard matrices with real diagonal of order 
6. Furthermore, relaxing the condition on the diagonal entries we obtain a new non-affine 
one-parameter orbit connecting the Fourier matrix F§ and Diva's matrix Dq. This answers 
a recent question of Bengtsson & al. in [3]. 
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1. Introduction 

Throughout this paper we will use the notation of [121 [T3] for well-known complex Hadamard 
matrices such as F 6 , D 6 , B 6 etc. 

Original interest in complex Hadamard matrices arose in connection with orthogonal pairs 
of maximal Abelian *-subalgebras (MASA's) of the n x n matrices [H [TOl [6j [H [9] . Subse- 
quently, it was realized in [14] that complex Hadamard matrices also play an essential role 
in constructions of teleportation and dense coding schemes in quantum information theory. 
This fact has given a new boost to the study of complex Hadamard matrices in recent years. 
On the one hand, several new and general constructions of such matrices have appeared 
IH EE21 El EEE] ■ On the other hand it is natural to try to fully classify complex Hadamard 
matrices of small order, as such characterization is currently available only up to order 5 
in [5]. Recently some progress has been made in the 6x6 case in [2] where all self-adjoint 
complex Hadamard matrices are characterized, and in [3] where numerical evidence is given 
of the existence of a conjectured 4-parameter family. While an algebraic form of such a 
4-parameter family (if it exists at all) remains out of reach, in this paper we present a previ- 
ously unknown non-affine one-parameter family of 6x6 complex Hadamard matrices which 
connects the Fourier matrix F e and Diva's matrix D 6 . This result complements the recent 
catalogue [T2] and answers a question of [3] , proving that apart from the isolated matrix S e 
the set of known 6x6 Hadamard matrices is connected. 

It is also important to mention that the 6x6 case is distinguished as the smallest dimen- 
sion where the maximum number of mutually unbiased bases (MUBs) is not known. It is 
well-known that if d is a prime power than the maximal number of MUBs in C d is d + 1. 
The existence of MUBs is equivalent to the existence of d x d complex Hadamard matrices 
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satisfying certain conditions (see e.g. [12] ). For the current status of MUB-related problems 
and, in particular, the case d = 6 we refer to [3] and references therein. The recent discovery 
of the new family in [2J, and the results of this paper may well be useful in the resolution 
of the MUB problem in dimension 6, and may give an indication to the maximal number of 
MUB's in dimension d = pq. 

Throughout the paper we restrict attention to dephased, symmetric complex Hadamard 
matrices of order 6 (the standard terminology dephased meaning the normalization condition 
that all entries of the first row and column are +1). It is quite natural to study the symmetric 
case for two reasons. First, the results of [2 J show that it is hopeful to obtain closed algebraic 
expressions if we require the matrix to satisfy certain symmetry assumptions (in [2J the 
self-adjoint case was classified). Second, the inspection of known 6x6 complex Hadamard 
matrices shows that many of them, such as F 6 , D 6 , C 6 , Sg, are equivalent to a symmetric one 
(throughout the paper we use the standard notion of equivalence (see e.g. [12]), i.e. H 1 and 
H 2 are equivalent, Hi = H 2 , if Hi = D1P1H2P2D2 with unitary diagonal matrices Di, D 2 and 
permutation matrices Pi,P 2 ). These two facts suggested that the set of symmetric Hadamard 
matrices of order 6 is on the one hand 'small' enough to be described in algebraic form and, 
on the other hand, 'rich' enough to contain interesting families of matrices. However, this 
intuitive approach turned out to be a little too optimistic in the first respect, and we needed 
to put further restrictions on the diagonal elements so that our algebraic calculations come 
to a comprehensible end. Accordingly, the outline of the paper is as follows. In Section [2] we 
fully classify dephased, symmetric complex Hadamard matrices of order 6 with real diagonal. 
It turns out that under this restriction well-known matrices emerge only. Therefore, in 
Section [3] we relax the condition on some diagonal entries and this leads to the discovery of 
a new non-affine one-parameter family. 

2. Symmetric matrices with real diagonal 

First we recall a simple but extremely useful result of [5] (see also [2J Lemma 2.6]). 

Lemma 2.1. Let u,v,s,t be complex numbers on the unit circle. Then 

(u + v)(s + t)(us + vt) el. □ 

We will also need the following elementary facts. In a dephased Hadamard matrix the sum 
of the entries in each row is (except for the first row where the sum is, of course, n). Given 
a row (xi, x 2 , x 3 , x 4 , x 5 , x 6 ) we distinguish two possibilities. First, if £ := x ^+ x ^+ x 3+ x i = q 
then X5 = — Xq. Second, and more importantly, if E 7^ and |S| < 1 then the last two 
coordinates are determined (up to change of order) as 

(1) x 5 ,6 = -S±i^| V / l-|S| 2 

The point is that — 2S = X5 + xq and it is easy to see geometrically that £5 and xq, being 
unit vectors, are determined as above. 

The main result of this section is the following 

Theorem 2.2. Let H be a dephased, symmetric complex Hadamard matrix of order 6 with 
real diagonal. Then H is equivalent to Sq or Dq . 

The proof is based on Lemma [2. 1[ and some considerations similar to those in [2J. 
The diagonal elements of H belong to { — 1, 1} by assumption. It is clear that there 
are either at least four l's in the diagonal or at most three. Therefore, after a possible 
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permutation of the rows and columns it is enough to consider the following two possibilities 
for the diagonal of H: 

(2) Diag(H) G {(1, 1, 1, 1, *, *), (1, -1, -1, -1, *, *)}, 

where the *'s stand for ±1. 

Lemma 2.3. Let H be a 6 x 6 symmetric complex Hadamard matrix of the form: 

111111 

1 1 x y * * 



(3) 



Then 



H 



1 x 1 z * * 

1 y z 1 * * 

1 * * * * * 

1 * * * * * 



(a) two ofx,y,z must be equal. 

(b) H is equivalent to Sq. 

(Note that we do not assume here that the last two diagonal entries are real; it is already 
implied by the above form.) 

Proof. First we prove (a). Let us denote = u,h 2 £ = f , ^3,5 = s, ^3,6 = t. We will use 
Haagerup's idea as in Lemma [2.11 By the orthogonality relations of rows 1, 2, 3 we have 



(4) 
(5) 



(6) 



2 + x + y = — (u + v) 
2+x+z= -(s + t) 
l + x + x + zy = — (su + tv). 



Now, Lemma [2.11 implies 



(7) (2 + x + y)(2 + x + z)(l + x + x + zy) G R. 

By similar arguments we obtain 



(9) 



(2 + x + y)(2 + y + z)(l + y + y + zx) G 
(2 + x + z){2 + y + z){l + z + z + yx) G 1 



After summing up these three expressions and eliminating real terms we get 



(10) 



x 2 y + xy z + xz 2 + x 2 z + y 2 z + yz 2 + 8 (xy + xz + yz) G 



Since a complex number is real if and only if it is equal to its conjugate, we obtain 
an equality if we replace each variable by its conjugate (i.e. its reciprocal) in the above 
expression. The resulting equality can then be rearranged by simple algebra to yield 



(11) (x — y) (x — z) (y — z) (xy + yz + zx + 8xyz + x 2 yz + xy 2 z + xyz 2 ^ = 



4 MATE MATOLCSI, FERENC SZOLLOSI 

The last factor in the product is clearly non-zero by the triangle inequality (one term has 
modulus 8, and the others have modulus 1). This proves (a). 



Now we turn to (b). It is easy to see that all x = y, 



x 



z, y = z lead to equivalent 



Hadamard matrices by permutation, thus we can assume without loss of generality that 
x = y. 

Now substitute back to (J7J) and eliminate all real terms to get 



(12) 



x + xz + xz G 



Therefore this expression equals its conjugate and we obtain 
(13) 



x 2 + xz + xz — x 2 — xz — xz = 



which yields 



(14) 



[x 2 — l) (z + x 2 z — x — xz 2 ") = 



Here x = y = 1 is clearly a contradiction, because the first two rows of H cannot be 
orthogonal. 

To show that x — y — — 1 is also impossible we need to consider two subcases. If 
x = y = — 1 = z then the rows of the leading 4x4 minor of H are already mutually 
orthogonal, therefore the last two entries of the first four rows of H should also be mutually 
orthogonal, which is clearly impossible, li x = y = —1 ^ z then H must be equivalent to a 
matrix of the following form 



(15) 



H 



- 1 


1 


1 


1 


1 


1 


1 


1 


-1 


-1 


u 


—u 


1 


-1 


1 


z 


—z 


-1 


1 


-1 


z 


1 


-1 


—z 


1 


u 


—z 


-1 


* 


* 


1 


— u 


-1 




* 


* 



The last two entries of rows 3, 4 are determined by (PQ), and the order of —1 and z in the 
fourth row is forced by the orthogonality of rows 3-4. The same orthogonality now implies 
z = io or z = uj 2 (with u being the third root of unity). If z = to then the orthogonality of 

implies u = again 



rows 2 and 3 implies that u = — A= which is a contradiction, z = cu' 
a contradiction. 

Therefore in (1T4T) we must have (z + x 2 z — x — xz 2 ) = (z — x)(l — xz) = 0, which implies 
z = x or z = x. The case z = x = y leads again to contradiction due to the following reasons. 
First, we cannot have x = y = z = — 1 as argued already above. Second, if x = y = z^— 1 
then (TjQ) implies that two of the rows 2, 3, 4 must contain the same two entries in the last 
two places. However, those two rows then contain the same element in four places, therefore 
they cannot be orthogonal. 
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Thus, we must have y = x ^ ±1, z = x. In this case H must be equivalent with the 
following form 

11111" 

1 x x * * 

X 1 X U V 
X X 1 V u 

* U V * * 

* V u * * 

where dt[x] < due to the orthogonality of rows 1 and 3. After conjugating H if necessary 
(and noting that Sq is equivalent to its conjugate) we can also assume that $S[x] > 0. The 
third and fourth rows are being forced by (JT]) as follows (we are free to choose the order due 
to permutation equivalence): 

(17) u = — 1 — U[x] + V-9?[a;] 2 -29?[x], 



(16) 



H 



(18) v = -1 - - V-SfW -2»[x]. 

Now, due to the orthogonality of rows 3, 4 we get 43ft [x] 2 + 109ft fx] +4 = which gives the 
only possible solution $t[x] = — \ and, by $$[x] > 0, we get x — u. Then by (iT7j) and ( TT8T) 

u — LO'V — u 2 , and 



(19) 



H 



1 
1 

uJ 



1 
1 



C<J id 



UJ 
2 



to 
u 2 
1 



1 

* 
UJ 
CO 2 

* 
* 



UJ 
UJ 

* 
* 



All the remaining entries are determined uniquely (using ([T]) and orthogonality) and we 
obtain 



(20) 



H 



1 
1 

UJ 



1 

UJ 

1 



1 1 

UJ UJ 2 



UJ UJ 



UJ 



UJ 



UJ UJ 



UJ 

1 

UJ 2 
UJ 



UJ 
UJ 2 

1 

UJ 



UJ 

UJ 2 
UJ 
UJ 

1 



This matrix is clearly equivalent to Tao's matrix Sq (we note here, that this matrix was 
published earlier in pQ, page 104). □ 



Having classified dephased, symmetric Hadamard matrices with diagonal (1,1,1,1,*,*) 
we can assume that the number of l's in the diagonal are at most three, in other words, the 
diagonal is (1, —1, —1, —1, *, *). 
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Lemma 2.4. Let H be a 6 x 6 symmetric complex Hadamard matrix of the form 



11 1 1 11" 

1 —1 x y * * 
1 x —1 z * * 
1 y z —1 * * 
1 * * * * * 
1 * * * * * 

Then H is equivalent to D e . (Note that we do not assume here that the last two diagonal 
entries are real.) 

Proof. We can assume that ^s[x] > (by conjugating all entries of H if necessary, and noting 
that Dq is equivalent to its conjugate). 

First assume that two of x, y, z are equal, say x = y (we are free to choose due to permu- 
tation equivalence) and, furthermore, z = —x. Then H is equivalent to 



(21) 



H 



- 1 


1 


1 


1 


1 


1 


1 


-1 


X 


X 


—x 


— X 


1 


X 


-1 


—x 


u 


—u 


1 


X 


—x 


-1 


V 


—v 


1 


—x 


u 


V 


* 


* 


1 


—x 


—u 


—v 


* 


* 



By the orthogonality of rows 3, 4 



(23) l + l + x + x + uv + uv = 0, 

which yields uv G R, therefore v — u or v — —u. The former is not possible due to (123]) . 
therefore v = —u. Then, using again ff23l) we get x + x = and hence, by the assumption 
5s[x] > 0, x — i. 

Then, using ([1]) we obtain that the last two entries of rows 5, 6 must be —1 and i in some 
order. The same order in both rows is not possible because then these rows would agree in 
four entries and could not be orthogonal. Therefore we have two possibilities: 





' 1 


1 


1 


1 


1 


1 ■ 




- 1 


1 


1 


l 


1 


1 " 




1 


-1 


i 


i 


— i 


— i 




1 


-1 


i 


i 


— i 


— i 


H = 


1 


i 


-1 


— i 


u 


—u 


, or H — 


1 


i 


-1 


— i 


u 


—u 








-1 












-1 








1 


i 


— i 


—u 


u 




1 


i 


— i 


— u 


u 




1 


— i 


u 


—u 


-1 


i 




1 


— i 


II: 


— it 


i 


-1 




1 


— i 


—u 


u 


i 


-1 




1 


— i 


—u 


u 


-1 


i 



In the first case, by the orthogonality of rows 4, 5 we get u = ±i. Both choices lead to 
Hadamard matrices equivalent to Dq. 

In the second case the orthogonality of rows 4, 5 yields u = ±1. Both choices lead to 
Hadamard matrices equivalent to D§. 

Let us now turn to the case when x = y, but z ^ —x. We will show that it is not possible. 
The last two entries of row 3 must be — x and — z and we are free to choose the order due to 
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permutation equivalence. Therefore, 



■ 1 


1 


1 


1 


1 


1 ■ 


1 


-1 


X 


X 


—x 


—x 


1 


X 


-1 


z 


—x 


—z 


1 


X 


z 


-1 


* 


* 


1 


—x 


—x 


* 


* 


* 


1 


—x 


— z 


* 


* 


* 



By the orthogonality of rows 2, 3 it we get 
(26) 1 — x — x + xz + 1 + xz = 0, 

thus x~z G R implying x = z or x = —z. The first case is not possible by ( |26l) . while the 
second contradicts our current assumptions. 

Now we turn to the case when x, y, z are all distinct. Again, we will show that it is not 
possible. We need to distinguish two subcases. First we assume that one variable is the 
negative of another, say y = —x (we are free to choose due to permutation equivalence). By 
assumption we cannot have z = x or z = —x (in the latter case z = y would hold). Therefore 
we have two choices to fill up the last two entries of rows 3, 4: 



(27) 





- 1 


1 


1 


1 


1 


1 ■ 




- 1 


1 


1 


1 


1 


1 ■ 




1 


-1 


X 


—x 


u 


—u 




1 


-1 


X 


—x 


u 


— It 


H = 


1 


X 


-1 


z 


—x 


—z 


or H — 


1 


X 


-1 


z 


—x 


— z 


1 












1 






-1 








—x 


z 


-1 


X 


—z 




—x 


z 


—z 


X 




1 


u 


—x 


X 


* 


* 




1 


u 


—x 


X 


* 


* 




1 


— u 


—z 


—z 


* 


* 




1 


—u 


—z 


—z 


* 


* 



In the first case the orthogonality of rows 3, 4 implies z — ±i and we can assume (by 
conjugation if necessary) that z — i. Then the orthogonality of rows 2, 3 and 2, 4 yield 
the equalities 1 — x — x — xi — ux + ui = and 1 + x + xi + x + ux + u\ = which, after 
summation, imply u = — i. Substituting back u = — i we get x = ±i. But this implies z = x 
or z = y(= —x) which contradicts our current assumptions. 

In the second case of ( 1271) the orthogonality of rows 3, 4 implies 1 — 1— ~z — z + x~z — xz = 0, 
which yields x~z — xz e R. This is only possible if x~z G R, that is x = z or — x = z = y, 
which are both excluded by assumption. 

Lastly, assume that all x, y, z are all distinct and none of them is the negative of another. 
By formula ([1]) there are the following four different possibilities to fill out the last two 
entries of rows 2, 3, 4 (we are free to fix the order of — x and — y in row 2 due to permutation 
equivalence): 



(28) 



1 

—x 
—x 

-y 

* 
* 



-y 

-z 
-z 



* 
* 



-x 
-x 

—z 
* 
* 



1 

-y 

—z 

-y 

* 
* 



-x 
—z 

-y 
* 

* 



-y 

— z 
* 
* 



-x 
—z 

—z 
* 
* 



-y 

-x 

-y 

* 
* 



Now we analyze these cases separately. 



8 MATE MATOLCSI, FERENC SZOLLOSI 

CASE[1]: consider the first possibility listed in (|28|) . that is 



(29) 



H 



- 1 


1 


l 


1 


1 


l ■ 


1 


-1 


X 


V 


—x 


-y 


1 


X 


-1 


z 


—x 


—z 


1 


y 


z 


-1 


-y 


—z 


1 




—x 


-y 


* 


* 


. 1 


-y 


—z 


—z 


* 


* 



Taking the inner product of rows 2, 3 

(30) 1 — x — x + yz + l+yz = 

Thus, yz G K and hence z = y or z = —y which are both excluded by assumption. 

CASE [2]: the second possibility listed in (1281) . Consider rows 2, 4 and apply the same 
simple argument as in CASE[1] above. 

CASE[3]: the third possibility listed in (|28|) . The orthogonality of rows 2, 3 implies 

(31) 1 — x — x + yz + xz + xy = 0, 

which means that xy + x~z + yz G R. Therefore this expression equals its conjugate, i.e. 

(32) xy + xz + yz = xy + xz + yz. 

Using that conjugates are the same as reciprocals this equation is equivalent to 

(33) (x + y)(x + z)(y-z) = 0, 

which contradicts our assumptions. 

CASE[4]: the fourth possibility listed in ( J28l) . Consider rows 3, 4 and apply the same 
simple argument as in CASE[1] above. □ 



3. A NEW FAMILY OF 6 X 6 COMPLEX HADAMARD MATRICES 

In Theorem 12.21 we have classified all dephased symmetric Hadamard matrices with real 
diagonal. Furthermore, from Lemmas 12.31 and 12.41 we see that it was, in fact, enough to 
specify four real entries in the diagonal. It is then natural to investigate the two remaining 
real options for the first four entries of the diagonal, i.e. the cases DiagH G (1, —1, 1, 1, *, *) 
and DiagH G (1, —1, —1, 1, *, *). 

Along the lines of Lemmas 12.31 and 12.41 a case-by-case argument shows that there exists no 
dephased symmetric complex Hadamard matrix with diagonal (1, —1, —1, 1, *, *). We do not 
include the details of this fruitless calculation. 

The last remaining case, DiagH G (1, —1, 1, 1, *, *), turns out to be the most interesting 
one. Unfortunately we do not have a full classification in this case, but we are able to obtain 
new matrices nevertheless. For some preliminary calculations we disregard the last +1 entry 
in the diagonal, and assume only that DiagH G (1, —1, 1, *, *, *). Due to the presence of the 
— 1 in the second row, the remaining entries must be x, —x and y, —y, and H takes the form 
(up to permutation equivalence) 
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(34) 



H 



- 1 


1 


1 


1 


1 


1 


1 


-1 


X 


-y 


V 


— X 


1 


X 


1 


a 


b 


c 


1 


-y 


a 


* 


* 


* 


1 


y 


b 


* 


* 


* 


1 


—x 


c 


* 


* 


* 



Using the orthogonality of rows 1, 2, 3 we apply Lemma [2. II as follows. 



(35) 
(36) 



x + y = -(-y-x), 
2 + x + b = - (a + c) 
x + x + by = — (—ay - 



cx) 



(37) 
therefore 

(38) (x + y) (2 + x + b) (1 -x + x + by) e R 
After expanding and eliminating the real entries one gets: 

(39) — 2x 2 — 2xy + 2xy + x 2 y — x 2 b — xyb + xyb + y + b E E 
This expression therefore equals its conjugate and by simple algebra we get 

(40) (y + x) [b 2 (l + x 2 ) + 6(2 - x 3 - 2x 3 y - 2x 2 - x + y + x 2 y + 2xy) - xy(l + x 2 )] = 

Unfortunately, we do not know how to handle the case when the second factor equals zero, 
therefore we need to settle for the simplifying assumption y = —x. We remark here, however, 
that there exist non-trivial solutions with y ^ —x too, such as the permuted version of 
Bjorck's cyclic matrix 



(41) 



1 

-1 
d 2 
-d 
d 
-d 2 



1 

d 2 
1 

-d 3 

-d 

d 2 



1 

-d 
-d 3 
-d 3 
-d 

d A 



1 

d 

-d 
-d 
d 

-1 



1 

-d 2 

d 2 

d 4 
-1 
-rf 4 



where d 



2 ' 1 V 2 • 

Having made the assumption y = —x the matrix H will now be determined up to permu- 
tation equivalence and possible conjugation. Now, H takes the form 



(42) 



H 



- 1 


1 


1 


1 


1 


1 ■ 


1 


-1 


X 


X 


—x 


—x 


1 


X 


1 


a 


b 


c 


1 


X 


a 


t 


u 


V 


1 


—x 


b 


u 


p 


Q 


. 1 


—x 


c 


V 


Q 


r 



From the orthogonality of rows 2, 3 and 1, 3 we get: 
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(43) 1 — x + x + xa — xb — xc = 

(44) l + x+l + a + b + c = 
Multiplying (14411 by x and then summing up and using x = 1/x we get 

,a*\ x 2 -2x-3 

(45) a = . 

This equation does have two solutions such that both x and a are on the unit ball, 



(Ap\ 1-V13 V-5 + 2V13 
(46) x lt2 = ± i 

Let us take x = x± (the other choice leads to the conjugate matrix), and hence 

7 -Vis .Vi9 + uVTs 

47 a = i 

K J 9 9 

Now, since 2 + x + a^0we can apply and obtain (up to change of order, which we are 
free to choose due to permutation equivalence) 



, N , -14 + 2^13 - V-58 + 34VT3 . V 134 + 22x/l3 - 8^-2446 + 730^ 
(48) b = 1- 



18 18 



. x -14 + 2VT3 + V-58 + 3AVTS . V134 + 22^ + 8^-2446 + 730^ 

(49) c— hi— 

v ; 18 18 

Next we find t, u and v. The orthogonality of rows 1, 4 and 2, 4 yield 

(50) 1 + x + a + t + u + v = 0, 

(51) 1 — x + ax + fx — ux — vx = 0. 
Multiplying (IBTl) by x and then summing up we obtain 

, , x 2 — 2x — 1 — 2a 

(52) t = . 

Substituting the values of x and a we get t — 1. 

Then, using ([1]) we obtain (the order being determined by the orthogonality of rows 3, 4) 

(53) u = c 

(54) v = b 

Finally we can use (OQ) once again to complete rows 5 and 6 as (the order being determined 
by orthogonality of rows 4, 5, 6): 



(55) p = r = 3- Vl3-iJ -21 + eVlS 



-19 + 4^ .2^-122 + 38^13 
(56) q = h 1 
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Therefore, we have obtained 



(57) M 6 



where x,a,b,c,p,q are determined by (H61) . (H7|) . (H8"j) . (T491 . (153]) and (I5"B"1) . respectively. 
It is easy to check that M 6 is indeed Hadamard. What we have shown above is that 
up to permutation equivalence Mq and its conjugate Mg are the only dephased symmetric 
Hadamard matrices with diagonal (1, —1, 1, *, *, *) and second row consisting of the elements 
X 2 1 f Ou . Ou - cc . 3j ^ . 

We will now proceed to show that M 6 is not contained in any of the previously known 
6x6 families. We need the following trivial 

Lemma 3.1. // a symmetric complex Hadamard matrix H is equivalent to a self-adjoint 
one, then it is also equivalent to its own conjugate i.e. H = H. 

Proof. Let if be a symmetric complex Hadamard matrix and suppose that it is equivalent to 
a self-adjoint one, say to A = A*. Then there are unitary diagonal Di, D 2 and permutational 
matrices Pi , P 2 such that 

(58) P 1 D 1 HD 2 P 2 = A = A* = P*D*H*DlP? 

By multiplying both sides with D 2 P 2 from the left and P\D\ from the right we get: 

(59) D 2 P 2 P X D X HD 2 P 2 P X D X = H* = H 

This clearly says that H = H. □ 

As a consequence we have 

Proposition 3.2. Mq and Mg are not equivalent to any previously known complex Hadamard 
matrix of order 6. 

Proof. We will use the Haagerup A-set of a matrix H = [hjk] , defined as 

(60) A H = {hijhkjhkihu for all 1 < i,j,k,l < 6}. 

It is well-known that Ah is invariant under equivalence (see [5]). 

The list of known 6x6 Hadamard matrices is as follows. The Fourier family (a, b) 

and its transposed ^Fq 2 \o, b)^j , the Dita family D^\c), and Tao's matrix are listed in 

[12] . The recently discovered non-affine family B^\a) is given in [2]. 

Mg is clearly inequivalent to Sq due to the Haagerup A-set being different. M 6 is inequiva- 
lent to any matrix in F 6 (a, b) since the third root of unity 10 G h. p {2), for every matrix in 
that family (i.e. for every a,b), while uj ^ A^/ 6 . The same is true for the transposed family 

^Fg 2 ^, b) j . M 6 is inequivalent to any matrix in D^\c) since i G ^~ D W^ f° r every matrix 
in that family, while i ^ A m 6 ■ Finally, using Lemma 13.11 M 6 is inequivalent to any of the 
matrices contained in the self-adjoint non-affine family B^\a) since Mq is inequivalent to 
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Mq, which can be seen as follows. Let a be as in M 6 . One can easily check that a 2 G Am 6 
and a 2 £ A M& , while a 2 Ajj 6 and a 2 e A^g . 

The same proof works for Mg, too. □ 

We will now generalize our discrete result above and construct a continuous one-parameter 
family stemming from M 6 . As M 6 is in block-matrix form it is quite natural to try and replace 
the l's on the main diagonal by some parameter d and consider dephased symmetric matrices 
of the following form: 

111 1 1 - 

1 X x —x —x 
x d a b c 
x a d c b ' 
x b c p q 
x c b q p 

The orthogonality of rows 2, 3 and 1, 3 imply 



(61) 



M 6 (x) 



(62) 1 — x + dx + ax — bx — cx = 

(63) l + x + d + a + b + c = 
Multiplying (162!) by x and then summing up we get 

(64) a = x 2 - 2 *- 1 _ d _ 

It is easy to see that < \x 2 — 2x — 1\ < A for each choice of x on the unit circle, therefore 
we can use to apply ([T]) to obtain the values of a and d as follows (we are free to choose 
the order due to permutation equivalence): 



T 2_9r-1 (x 2 -2x-l) a/16 - \x 2 - 2x - 1| 2 

(65) a = ^— ^ - i- ^- — - 

v ; 4 4 |x 2 — 2x - 1| 

x 2 _9x-l (x 2 - 2x - 1) \/l6- |x 2 -2x- II 2 

(66) d = + i- ^4 l — - 

K ' 4 4|x 2 -2x-l| 

From (jMj) we see that 1 + x + a + d = which vanishes if and only if x = ±i. We exclude 
x = ±i from these considerations and remark that this case can be handled separately as in 
Lemma \2 .41 and leads to H = D 6 . For s ^ ±i the values of b and c are determined uniquely 
by ([I]) as follows (we are free to choose the order due to permutation equivalence): 



1+X 2 _(1 + x 2 )V16-|1+^ 
(67) b = i 



411 + x 2 \ 



(68) c = — + i 



l +x 2 ^1 + X *)JIQ-\1+ X 2\ 2 



4 4 11 + x 2 \ 
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It is easy to check (rather by computer) that with these parametric choices the first four 
rows of H are mutually orthogonal to each other. 

We evaluate 1 — x + 6 + cin order to use ([T]) again to determine p and q. The orthogonality 
of rows 1, 5 and 2, 5 imply 

(69) l-x + b + c + p + q = 

(70) 1 + x + xb + xc — xp — xq = 0. 
Multiplying (1701) by x and summing up we get 

-x 2 - 2x + 1 



(71) 



1 — x + b + c 



and we see that it does not vanish for unit vectors x. Therefore we can apply ([I]) to determine 
p and q as follows (the order now being forced by orthogonality): 



(72) 



(73) 



V 



x i + 2x-l ,( x + 2x - 1) W16 - \x 2 + 2x-l 



4 A\x 2 + 2x - 1| 

x 2 + 2x - l _ (x 2 + 2x - 1) yjl6- \x 2 + 2x - 1| 



— i 



4 " 4\x 2 + 2x-l\ 

It is easy to check by computer that the emerging matrix H is Hadamard for all i ^ ±i. 
Therefore we have obtained the following 

Theorem 3.3. There is a one parameter symmetric non-affine family of Hadamard matrices 
given by (ISTj) . with x = e**, x ^ ±i, and a, b, c, d,p, q being given as in (165|) . (167|) . (168|) . (1661) . 
(172|) . (1731) . respectively. 

Finally we make the following observation which answers a question raised in [3] and shows 
that the set of currently known Hadamard matrices of order 6 is connected except for the 
isolated matrix 5*6- In particular, the family M 6 (x) above connects the Fourier matrix F 6 
and Dita's matrix D§. 



Observation 3.4. M 6 (l) = F 6 , lim M 6 (e a ) = D 6 , and finally, M 6 ( e 



(74) 



(75) 



lim M 6 (e 1 *) 
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In summary, we have given a full classification of dephased symmetric complex Hadamard 
matrices of order 6 with real diagonal, showing that in this case the well-known matrices Sq 
and Dq emerge only Furthermore, relaxing the reality condition on the diagonal entries we 
have been able to obtain a new non-affme family of dephased symmetric Hadamard matrices 
of order 6 which connects Fq and D§. 

It would be interesting to see whether this family can be extended by further parameters. 
For example, it is natural to try to replace the second entry of the diagonal (—1) by a 
parameter h. However, currently we are unable to classify that case. It also remains to be 
seen whether this new family helps to increase the number of bases appearing in MUBs in 
C 6 . 
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